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Abstract.  Films  and  molds  of  nematic  polymer  materials  are  notorious  for  heterogeneity  in  the 
orientational  distribution  of  the  rigid  rod  or  platelet  macromolecules.  Predictive  tools  for  structure 
length  scales  generated  by  shear-dominated  processing  are  vitally  important:  both  during  processing 
because  of  flow  feedback  phenomena  such  as  shear  thinning  or  thickening,  and  postprocessing  since 
gradients  in  the  rod  or  platelet  ensemble  translate  to  nonuniform  composite  properties  and  to  residual 
stresses  in  the  material.  These  issues  motivate  our  analysis  of  two  prototypes  for  planar  shear 
processing:  drag-driven  Couette  and  pressure-driven  Poiseuille  flows.  Hydrodynamic  theories  for  high 
aspect  ratio  rod  and  platelet  macromolecules  in  viscous  solvents  are  well  developed,  which  we  apply 
in  this  paper  to  model  the  coupling  between  short-range  excluded  volume  interactions,  anisotropic 
distortional  elasticity  (unequal  elasticity  constants),  wall  anchoring  conditions,  and  hydrodynamics. 
The  goal  of  this  paper  is  to  generalize  scaling  properties  of  steady  flow  molecular  structures  in  slow 
Couette  flows  with  equal  elasticity  constants  [M.  G.  Forest  et  al.,  J.  Rheol. ,  48  (2004),  pp.  175—192] 
in  several  ways:  to  contrast  isotropic  and  anisotropic  elasticity;  to  compare  Couette  versus  Poiseuille 
flow;  and  to  consider  dynamics  and  stability  of  these  steady  states  within  the  asymptotic  model 
equations. 
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tions,  instability 
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1.  Introduction.  Shear  dominated  flows  of  nematic  liquid  crystal  polymers 
(NLCPs)  generate  anisotropy  and  spatial  heterogeneity  in  the  orientational  distri¬ 
bution  of  the  rigid  rod  or  platelet  ensemble.  These  phenomena  are  well  documented 
in  light  scattering  textures  [9,  1,  24,  25,  31].  A  characterization  of  the  lengthscales  in 
the  molecular  distribution  responsible  for  the  scattering  patterns,  and  whether  they 
are  due  to  changes  in  the  direction  of  peak  orientation  (nematic  elasticity)  or  due  to 
focusing  and  defocusing  of  the  orientational  distribution  (molecular  elasticity),  are 
the  subject  of  numerous  modeling  and  computational  studies  (cf.  [32,  30,  21,  29]). 
Molecular  orientation  features  in  different  flow  regimes  are  of  extreme  importance 
for  materials  design,  as  they  impart  anisotropic  and  nonuniform  material  properties 
[34,  18,  19].  Another  issue  typical  of  non-Newtonian  fluids  is  flow  feedback,  where 
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elastic  stresses  alter  apparent  viscosity.  This  paper  is  a  continuation  of  our  systematic 
studies  of  mesostructures,  both  from  free  space  elasticity  patterns  absent  of  external 
fields  and  boundary  anchoring  conditions  [11,  12],  and  from  planar  Couette  cells  mov¬ 
ing  at  prescribed  slow  speeds  [13,  14].  We  defer  to  these  articles,  where  a  detailed 
account  of  analytical  results  on  continuum  Leslie-Ericksen-Frank  (LEF)  models  is 
given,  notably  by  [26,  4,  5,  7,  27,  23,  28]. 

In  [14],  the  authors  considered  a  Doi-Marrucci-Greco  (DMG)  mesoscopic  orienta¬ 
tion  tensor  model,  allowing  a  full  coupling  between  flow  structure,  director  (nematic) 
and  order  parameter  (molecular)  distortions,  and  with  imposed  plate  motion  and 
molecular  anchoring  conditions.  The  model  has  been  benchmarked  in  the  longwave, 
monodomain  regime  with  resolved  simulations  of  the  Doi  kinetic  theory  [15,  16].  In¬ 
deed,  the  motivation  for  an  analytical  study  of  structure  properties  is  to  provide  guid¬ 
ance  for  structure  simulations  of  mesoscopic  [17]  and  kinetic  [15,  16]  models,  where 
the  parameter  space  is  too  large  to  assimilate  any  kind  of  collapse  of  the  numerical 
data  through  scaling  laws. 

In  this  paper,  we  extend  our  previous  asymptotic  scaling  analysis  in  several  ways. 
First,  we  consider  a  more  general  physical  model  to  admit  anisotropic  distortional 
elasticity  (unequal  bend,  splay,  twist  elasticity  constants).  The  second-moment  orien¬ 
tation  model  is  derived  from  a  recent  generalization  of  the  Doi-Hess-Marrucci-Greco 
kinetic  theory  [33]  and  guides  our  numerical  studies  [15,  16],  for  which  there  are 
no  preceding  numerical  or  analytical  results.  Second,  the  asymptotic  analysis  is  ex¬ 
tended  from  plate-driven  Couette  cell  properties  to  pressure-driven  Poiseuille  flows. 
The  boundary  conditions  consist  of  molecular  orientational  anchoring  conditions  at 
solid  walls,  where  the  degree  of  order  is  set  by  the  concentration  of  the  nematic  liquid 
and  the  principal  orientation  axis  is  a  free  parameter,  together  with  no-slip  conditions 
for  the  velocity  field.  We  further  assume  an  in-plane  orientation  tensor  (restricting 
the  principal  orientation  axes  of  the  molecular  distribution  to  the  flow-flow  gradient 
plane),  and  posit  that  the  velocity  field  varies  only  transverse  to  the  primary  flow  di¬ 
rection.  These  assumptions  are  not  easily  lifted,  in  that  the  fortuitous  diagonalization 
of  the  flow-nematic  steady  balance  equations  is  apparently  lost  for  higher  dimensional 
orientational  and  spatial  degrees  of  freedom.  Finally,  we  extend  the  asymptotic  anal¬ 
ysis  to  time-dependent  model  equations. 

From  this  formulation,  we  develop  a  formal  asymptotic  analysis  in  the  slow-plate 
(so-called  small  Deborah  number)  and  weak  pressure  gradient  limits,  which  yield 
exactly  solvable,  steady  flow-nematic  model  equations.  From  the  explicit  solutions, 
lengtliscale  selection  criteria  and  scaling  properties  become  explicit,  parameterized 
in  terms  of  molecular  parameters  (nematic  concentration  N,  molecule  aspect  ratio  r, 
persistence  length  C  of  distortional  elasticity,  persistence  length  L  of  the  anisotropic 
distortional  elasticity)  [33],  and  experimental  conditions  (gap  width  (2 h),  plate  speeds 
±wo  for  plane  Couette  flows  and  pressure  gradient  for  the  plane  Poiseuille  flow, 
and  plate  anchoring  conditions  on  the  molecular  field).  From  the  time-dependence  in 
the  asymptotic  equations,  we  explore  transient  solutions  at  the  first  and  second  order 
in  the  asymptotic  scheme  to  infer  stability  of  the  steady  states  within  the  asymptotic 
balance  equations.  We  first  consider  plane  Couette  flow,  followed  in  the  next  section 
by  plane  Poiseuille  flow. 

2.  Spatial  structures  and  their  stability  in  plane  Couette  flows.  We 

consider  plane  Couette  flow  between  two  parallel  plates  located  at  y  =  ±h  and  moving 
with  velocity  v  =  (±uo,0,  0),  respectively,  in  Cartesian  coordinates  ( x,y,z ).  Figure  1 
depicts  the  cross  section  of  the  flow  geometry  on  the  ( x ,  y)  plane. 

Here  we  consider  flow-orientation  interactions  in  weak  plane  Couette  flow,  char- 
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v=-v  0  >  Q=Q0 

Fig.  1.  Geometry  of  plane  Couette  flow.  The  gap  width  in  the  shear  cell  is  2h.  The  LCP  in 
the  cell  is  sheared  by  moving  the  upper  plate  with  a  constant  speed  vq  and  the  lower  one  with  the 
same  speed  in  the  opposite  direction.  At  the  bounding  surfaces,  the  orientation  tensor  is  equal  to  its 
equilibrium  value. 


acterized  by  a  small  effective  or  averaged  shear  rate.  We  nondimensionalize  using 
the  gap  half-width  (h)  between  the  shearing  plates  and  the  nematic  polymer  mean 
relaxation  time  to  =  yjo  7  where  D®  is  the  rotary  diffusivity  for  the  rigid  rod  or  platelet 
[33].  We  denote  the  position  vector  by  x,  the  velocity  by  v,  the  extra  stress  tensor 
by  r,  and  the  pressure  by  p,  respectively.  The  dimensionless  flow  and  stress  variables 
are  defined  by: 


(1) 


_  h 2 


P  =  ~rP, 
Jo 


where  /o  =  phA  /£$  is  a  mesopliase  bulk  force  and  p  is  the  nematic  polymer  (NLCP) 
density.  Let  c  be  the  NLCP  number  density,  k  the  Boltzmann  constant,  T  absolute 
temperature,  N  a  dimensionless  concentration,  ps  the  solvent  viscosity,  and  Q,i  = 
1,  2,  3  three  friction  coefficients  related  to  NLCP-solvent  interactions.  £  measures 
the  range  of  isotropic  elastic  interaction  while  L  does  so  for  the  anisotropic  elastic 
interaction  [33].  The  following  eight  dimensionless  parameters  arise: 


(2)  Re 


ph2  _  3 ckTt0 

toys  ’  h2p  ’ 


Er  = 


8ft2 
NC2  > 


Pi  = 


3ckT£ito 
Wp  ’ 


*  =  1,2,3,  9  =  %. 


a  measures  the  strength  of  elastic  energy  relative  to  kinetic  energy;  Re  is  the  solvent 
Reynolds  number;  Er  is  the  Ericksen  number  which  measures  the  relative  strength  of 
the  short-range  nematic  potential  and  the  isotropic  distortional  elasticity  potential;  9 
measures  the  degree  of  anisotropy  in  the  distortional  elasticity,  with  values  limited  to 
[—1,  00);  1  / pi,  i  =  1,  2,  3  are  three  nematic  Reynolds  numbers.  We  drop  the  tilde'  on 
all  variables  from  now  on  so  that  all  equations  and  figures  in  the  following  correspond 
to  normalized  variables,  length,  and  time  scales. 

The  dimensionless  forms  of  the  balance  of  linear  momentum,  stress  constitutive 
equation,  and  the  continuity  equation  (dimensional  forms  are  in  [33] )  take  the  follow¬ 
ing  form. 

Linear  momentum  balance: 


(3) 


|v  =  V.(-pI  +  r), 


where  external  forces  are  neglected. 
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Continuity  equation: 


(4)  V  •  v  =  0. 

Constitutive  equation  for  the  extra  stress: 

r  =  2?;D  +  aa[M  -  §  -  f  ((I  +  ^  A)M  •  M  +  M  •  (I  +  ^A)M 
-2(1  +  ^A)M  :  M4)]  -  ^(AM  •  M  -  M  •  AM)  -  ^[VM  :  VM 

—  (VVM)  :  M]  +  f^[4M6  ::  VVM  +  2M4VV  ::  M4  -  VVM:M4  -  (VVM:M4)T 
— M4:VVM  -  (M4:VVM)t  -  (MVV:M4)T  -  MVV:M4]  -  Tf|F[VVM:M4 

—  (VVM:M4)t  -  M4:VVM  +  (M4:VVM)t  -  MVV:M4  +  (MVV:M4)T] 

+[fii(a)(DM  +  MD)  +  /X2(cj)D  :  ]VI4], 

(5) 

where  M  is  the  second  moment  of  the  orientational  probability  density  function  of 
the  kinetic  theory  (called  the  structure  tensor),  M4  and  Mg  are  the  fourth  and  sixth 

moment  of  the  probability  density  function,  respectively,  ?y  =  1/Re  +  ^3(0)  and 
2_|_1 

a  =  ■/tzj  parameterizes  the  aspect  ratio  r  of  the  spheroidal  molecules,  where  0  <  a  <  1 
corresponds  to  a  rod- like  molecule  and  —  1  <  a  <  0  for  platelets  [33]. 

The  boundary  conditions  on  velocity  v  are  scaled  to 

(6)  v|y=±4  =  (±De,0,0), 


where 


(7) 


the  Deborah  number ,  is  the  ratio  of  the  relaxation  time  relative  to  the  time  scale  set 
by  the  moving  plates  in  the  shear  experiment.  Weak  shear  is  defined  by  a  small  value 
of  De  indicating  the  time  scale  set  by  the  shear  experiment  is  much  larger  than  the 
molecular  relaxation  time  scale.  Following  previous  studies  [8,  11,  12],  we  assume 
strong  molecular  anchoring  at  the  plates  given  by  the  quiescent  nematic  equilibrium 
of  the  orientation  tensor  (the  deviatoric  part  of  the  structure  tensor)  Qo  =  Mo  —  |  = 
So(nn  —  |).  The  rest  state  equilibrium  of  Q  at  sufficiently  high  concentrations  is  a 
uniaxial  nematic  phase,  with  unique  order  parameter, 


(8) 


1 

/ 

8 

— 

1+3v 

1 - 

4 

3  N 

The  uniaxial  director  n  is  arbitrary  for  quiescent  phases;  this  degeneracy  is  bro¬ 
ken  experimentally  by  mechanical  or  chemical  plate  preparations.  We  model  a  uni¬ 
form  plate  anchoring  condition,  either  parallel  to  the  flow  direction,  called  tangential 
anchoring,  or  perpendicular  to  the  shearing  plates,  called  normal  (or  homeotropic) 
anchoring. 
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The  time  evolution  equation  of  M  (in  dimensionless  form)  is  given  by  [33] : 

±M  —  fi  •  M  +  M  •  fi  —  a[D  •  M  +  M  •  D]  =  -2aD  :  M4 
-6[Q  -  N(M  •  M  -  M  :  M4)]  +  ^[AM  •  M  +  M  •  AM  -  2AM  :  M4] 

+  2L  [(VVM)[M4  +  ((VVM)]M4)T  +  m4!vvm  +  (m4!vvm)t 

+MVV:M4  +  (MVV:M4)t  -  4M6  ::  VVM  -  2M4VV  ::  M4], 

In  order  to  arrive  at  a  closed  system  of  governing  equations  at  the  level  of  second 
order  tensors,  we  approximate  fourth  (M4)and  sixth  (Mg)  order  tensors  in  the  above 
governing  system  of  equations  using  the  following  simple  closure  rules: 

(10)  M4  « MM,  M6  « MMM. 

These  simple  closure  approximations  respect  the  traceless  property  of  the  orientational 
dynamic  equation,  and  have  been  shown  to  yield  a  good  approximation  of  kinetic 
theory  in  the  dynamics  of  monodomains  at  the  nematic  concentrations  of  interest 
here  [10,  15,  16].  These  closures  are  exact  when  the  molecules  are  aligned  perfectly. 

We  remark  that  the  distortional  elastic  free  energy  reduces  to  the  Oseen-Frank 
energy  after  the  closure  approximation,  in  which  the  three  Frank  elastic  constants  are 
given  by 

(11)  =  k2  =  l^s2(  1  +  f  (1  -  s)),  k3  =  f^2(l  +  |(1  +  4s)). 

For  rod-like  NLCPs, 

(12)  0  <  k\  =  k2  <  kz; 
whereas  for  discotic  NLCPs  (platelets), 

(13)  0  <  kz  <  k\  =  k2 . 

2.1.  Asymptotic  solutions  in  weak  plane  Couette  flows.  We  seek  asymp¬ 
totic  solutions  of  the  governing  system  of  equations  with  the  boundary  conditions 
given  by  (6)  and  (8).  We  employ  a  biaxial  representation  of  the  orientation  tensor 
[11] 

(14)  Q  =  s  |nn  -  ^lj  +  0  ^nxnx  -  ^iV 

where  (s,  (3)  are  two  order  parameters  measuring  the  birefringence  relative  to  the 
optical  axes  (also  called  directors)  n  and  nx  confined  to  the  shearing  plane  ( x ,  y)  and 
parameterized  by  a  director  angle  ip, 

(15)  n  =  (cosip,  sin  ip,  0),  nx  =  (—  sin?/’,  cos  ip,  0), 
and  I  is  the  3x3  identity  matrix.  We  propose  the  solution  ansatz 

OO  OO  OO 

(16)  vx  =  (•)  =  (9)kDek,  ip  =ipo  +  '^/i>{k)Dek, 

fc= 1  fc=0  fc= 1 

where  (•)  represents  the  order  parameters  s,  (3,  respectively.  The  solution  is  sensitive 
to  the  choice  of  boundary  conditions,  so  we  present  tangential  pipo  =  0)  and  normal 
(ipo  =  f )  anchoring  conditions  separately. 
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2.2.  Tangential  anchoring  (ipo  =  0).  First,  we  note  that 


(17)  v£k)  =  ^(2fc)  =  s2fc— 1  =  /?2fc— l  =  0,  k  =  1, . . . ,  00, 


demanded  by  the  boundary  conditions  and  the  governing  equations  at  the  respective 
orders.  This  also  applies  to  the  case  of  normal  anchoring,  but  not  to  tilted  anchoring 
(ipo  ^  0,  f )  [14].  The  governing  equations  at  order  0(1)  give  the  equilibrium  solution 
of  Q  consistent  with  the  boundary  anchoring  condition;  the  equations  at  order  O(De) 
are  obtained  by  solving  the  following  equations  for  ipA)  and  vi1'1: 


ch/jW  8 2ipd) 

dt  dy 2  dy 


(18) 


where 


(19) 


dvi1]  dr, 
dt 


‘xy 

dy 

_  92V’(1)  ,  n5vi1) 

Txv-C~d^+D^y' 


A=^.(s0  +  2)(3  +  8(l-s0)),  B  =  \(  Al-1), 


C=- 


18  Er 


[0{  1  -  S0)XL  +  3(Al  -  1)],  D  =  Um So  +  3 77), 


where  the  “tumbling  parameter”  A^  is  defined  by 

a(2  +  sq) 


(20) 


A  L  = 


3s0 


A/-J  >  1  corresponds  to  flow  aligning  and  |Al|  <  1  yields  director  tumbling  in  mono¬ 
domain  shear  flows  (Er  — >  00)  [14]. 

2.2.1.  Steady  state  features  of  the  major  director  •j/d1)  and  primary  flow 

The  nonzero  leading  order  steady  solution  for  the  velocity,  order  parameters  s 
and  (3 ,  and  the  director  angle  ip  can  be  solved  explicitly: 


(2!)  vi1)(y)  =  y,  s0  =  s0,  /?o  =  0,  ipw(y)  =  MEr(y2  -  1), 


(22)  M  4(s0+2)(3+(i-so)e) 

Note  that,  as  in  the  isotropic  elasticity  limit  [14],  these  solvability  conditions  imply 
simple  shear  flow  at  leading  order  in  De ,  and  yield  that  the  orientational  distribution 
is  dominated  by  nematic  (director)  distortions.  The  prefactor  (22)  yields  that  the 
winding  number  of  the  major  director  between  the  plates  is  proportional  to  the  Er- 
icksen  number,  as  with  the  isotropic  elasticity  limit  [14].  The  formula  (22)  yields  the 
scaling  law  for  elastic  distortions  which  are  nonuniform  across  the  gap  with  length- 
scale  proportional  to  M^1 ,  which  in  turn  is  proscribed  by  three  material  parameters: 
( a,N,6 ).  For  fixed  ( a,N ),  |M|  decreases  as  8  increases  (8  £  [— l,oo)).  We  summarize 
the  dependence  of  M  on  8  for  given  material  parameters  (a,  N)  as  follows: 

•  The  sign  of  M  governs  the  “chirality”  of  nematic  distortion,  or  direction  of 
director  rotation  from  the  plates.  M  is  negative  for  flow-aligning  rods  (a  > 
0,  A l  >  1)  and  positive  for  tumbling  rods  (a  >  0,0  <  Xl  <  1)  and  discs  or 
platelets  (a  <  0).  |M|  decreases  with  respect  to  all  8  £  [— l,oo). 
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Fig.  2.  M  as  a  function  of  N  and  6  with  tangentially  anchored  boundary  conditions  in  weak 
Couette  flow:  The  top  left  panel  is  for  rods,  with  a  =  0.8,  whereas  the  top  right  panel  is  for  platelets, 
with  a  =  —0.8.  The  bottom  panels  show  M  for  two  distinct  concentrations  of  rods  (left)  and  platelets 
(right). 


•  Physically,  the  director  angle  winds  counterclockwise  for  rods  in  the  flow¬ 
aligning  regime  Xl  >  1  from  the  lower  plate  to  the  midplane,  then  unwinds 
from  the  midplane  to  the  upper  plane;  the  orientation  reverses  for  rods  in  the 
tumbling  regime  and  for  platelets  in  all  regimes. 

•  At  fixed  Ericksen  number,  anisotropic  elasticity  tends  to  reduce  the  magnitude 
of  the  director  winding  for  rods  while  enhancing  director  distortion  for  platelets. 

Figure  2  depicts  M  as  a  function  of  ( N ,  6)  two  values  of  a  corresponding  to  rods 
(a  =  0.8)  and  platelets  (a  =  —0.8).  Next,  we  consider  a  limited  notion  of  stability  of 
this  steady-state  structure,  by  studying  transients  of  the  first  order  governing  system 
of  equations  in  the  presence  of  superimposed  spatial  disturbances. 

2.2.2.  Transient  behavior  of  (id1) ,  i/d1))  near  steady  states.  The  transient 
solution  for  v i1'*  and  i/d1)  (the  difference  between  the  time-dependent  solution  and  the 
steady  state)  obeys  the  same  homogeneous  linear  partial  differential  equations  but 
satisfies  a  zero  boundary  condition.  Its  behavior  dictates  the  stability  of  the  steady 
state  within  the  asymptotic  balance  model:  the  steady  state  is  asymptotically  stable 
if  the  transient  solution  vanishes  as  t  — >  oo. 

Proposition  1.  The  steady  solution  (21,  22)  of  system  (18)  is  stable  for  AD  — 
BC  >  0  and  unstable  for  AD  —  BC  <  0  with  respect  to  zero  boundary  conditions  on 
and  i/d1). 
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Proof.  We  first  consider  the  case  of  BC  ^  0  and  prove  the  steady  solution  is 
stable  provided  AD  —  BC  >  0.  In  the  following  proof,  we  drop  the  superscripts  on  if 
and  vx.  Extending  (18)i  to  the  boundary  and  accounting  for  the  boundary  condition 
if(— 1,  t)  =  tf(  1,  t)  =  0,  we  have 

(23)  (^  +  Bt)U=0. 

We  introduce  a  nonnegative  functional 

(24)  I(t)  =  J  [8 iifl  +  S2 vl\dy 

with  <5i  >  0  and  82  >  0.  We  note  that  A  >  0  and  D  >  0  from  (19).  □ 

Case  1.  BC  <  0.  Choosing  8\  =  \C\  and  82  =  \B\  and  integrating  by  parts,  the 
time  derivative  of  the  nonnegative  functional  can  be  estimated: 

^  ^  =  “2  J  [SiAifyy  +  (SiB  +  82C)ifyyvX}y  +  82Dv1y]dy 

(25) 

=  — 2  J  [\C\Aifyy  +  \B\Dv2xJdy  <  0. 


This  shows  that  the  steady  solution  of  the  system  is  stable. 

Case  2.  Choose  £1  =  max(^U,  1)  >1,62  =  max(|n,  1)  >  1.  We  have 


dl(t) 

dt 


=  2 


(26) 


1-1 


[SlifyAy  +  S2VxVxt\dy 


=  -2 


/-1 


[biAifl  +  ( 81B  +  b2C)ipyyvx,v  +  b2Dv2x  }dy. 


The  integrand  is  quadratic  and  the  discriminant  is 
{biB  +  b2C)2  -  AbxAb-zD 


(27) 


+  82AD 


2b1b2{AD  -  BC) 


=  biAD 


B 2  (  C2 

—r —  max  — — 

AD  \AD 


,11—  max  (  ,  1 


/  B2 
\  AD  ’ 


-\-b2AD 


C2 

AD 


max 


I  1 


-2bib2{AD  -  BC)  <  —2b\82{AD  -  BC)  <  0. 

t~>  2  f~i2 

The  first  inequality  is  based  on  <  1  because  of  AD  —  BC  >  0  and  BC  >  0. 

Since  b\A  >  0,  the  integrand  is  always  positive;  thus  <  0.  Hence  the  steady 
solution  of  the  system  is  stable. 
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The  proof  for  BC  =  0  is  far  simpler  and  omitted.  To  prove  instability  when 
AD  —  BC  <  0,  we  only  need  to  find  one  unstable  mode.  Let  )  =  ipy(y,t)',  the 
system  (18)  with  the  boundary  condition  becomes 

( fit  —  A(f)yy  +  Bvx,yyi 


(28)  Vx,t  —  ^ 4*yy  ~b 


(A(f)y  +  Bvx^y)\y~ -(- 1  0, 


vx(l,t)  =  1,  vx(-l,t)  =  -l. 


To  find  an  unstable  mode,  we  seek  normal  modes  of  the  form 


(29) 


ky)  ^ 
vx{y)  ) 


and  consider  the  resultant  eigenvalue  problem 


(30) 


P  = 


A  B  \ 
CD  )' 


where  7  is  the  growth  rate.  The  steady  solution  is  unstable  if  7  >  0. 
The  matrix  P  has  two  distinct  eigenvalues  given  by 


(31)ri  =  A+D-^/{A+D)2-4{AD-BC)  <  Q 
Let  £  =  arctan  '4^ri ,  then 


_  A+D+J(A+D)2-4(AD-BC) 

r2  -  - o -  >  u 


(32) 


cos  £  sin  £ 
—  sin  £  cos  £ 


\ 

/ 


is  a  solution  of  (28)  satisfying  vx(l,  t )  =  vx(—  1,  t)  =  0  and  7  is  determined  by  ( A<j)y  + 
BvXty)\y-±i  =  0,  which  yields 


(33)  {A  -  ri)(l  +  At  g2fl) )  coth  -  ^-nr2  cot  =  0. 


As  7  — >  +00  ,  the  first  term  goes  to  a  finite  value  while  the  second  one  is  periodic 
and  varies  between  —00  and  +00  within  one  period.  Consequently,  the  equation  has 
infinitely  many  positive  solutions  for  7,  which  completes  the  proof,  and  indicates  that 
the  diagnostic  AD  —  BC  signals  catastrophic  instability  when  it  is  negative. 

We  note  that  for  discs  (a  <  0)  and  flow-aligning  rods  (a  >  0,  \l  >  1), 

AD  —  BC  =  2Yet [3(2  +  soXatsq  +  3 rj)  +  f  ((AL  -  l)3s0)2 

(34) 

+0(1  —  so)(2  +  so)(^a((AL  —  l)3so)  +  piso  +  3ry)]  >  0. 

Hence,  the  steady  state  is  always  stable  for  discotic  LCPs  and  flow-aligning  rods,  and 
may  be  unstable  only  for  tumbling  rods  (0  <  a,  0  <  Xl  <  1). 
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Fig.  3.  The  neutral  stability  curves  defined  by  A1B2  —  A2B1  =  0  and  AD  —  BC  =  0  with 
tangentially  anchored  boundary  conditions  in  weak  shear.  The  arrows  point  to  the  positive  directions 
of  the  discriminant.  The  systems  for  (s2,^2)  and  are  stable  when  the  discriminant  is 

positive,  respectively,  unstable  otherwise.  The  parameter  values  are  Er  =  100,  a  =  1,/xi  =  0.001 
and  rj  =  0.002. 


The  neutral  stability  curve  {AD  —  BC  =  0)  depends  on  the  values  of  the  param¬ 
eters  (0,  a ,  N,  a ,  /xi,  77).  From  (34),  AD  —  BC  >  0  can  be  translated  into  two  separate 
constraints  on  the  energy  parameter  a  and  0, 


(35) 


a  <  ac 


4(ahsq+3?))  . 
off—  \l)so  ’ 


a  >  ac  and  9  <  9C 


3(2+so)(/iiso+3??)  +  ^  —  l)3so)2 

(l-so)(2+so)(Ho((l-AI,)3so)-7iiso+3r/)  ’ 


Figure  3  depicts  the  stability  transition  curve  in  the  parameter  space  (N,  9)  at  a 
few  selected  values  of  other  parameters.  We  observe  that  the  values  of  N,  9  that 
yield  instability  tend  to  be  large  and  out  of  the  practical  range  for  nematic  polymer 
materials.  We  also  note  that  the  flow-aligning  region  and  the  stable  region  versus 
(N,  9)  both  grow  significantly  as  the  shape  parameter  a  increases,  i.e.,  as  the  aspect 
ratio  becomes  more  extreme.  The  instability  region  vanishes  as  a  — »  The  stable 

region  also  grows  as  7 ii/a,ri/a  increase. 

We  summarize  this  more  precise  statement  of  Proposition  1  in  the  following  corol¬ 
lary. 

Corollary  1.  The  steady  asymptotic  solution  (21,  22)  is  asymptotically  stable 
within  the  leading  order  balance  equations  if  either  of  conditions  (35)  is  satisfied.  If 
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a  >  ac  and  9  >  9C,  where  aCl9c  are  defined  in  (35),  the  steady  solution  is  unstable 
and  the  leading  order  system  of  equations  is  ill-posed. 

2.2.3.  Steady  state  features  of  the  order  parameters  S2  and  /32  at  0(De2). 

The  order  parameters  (s,  /3)  vanish  at  leading  order  in  He,  with  order  0(De2)  behavior 
governed  by  the  equations 


(36) 


JH  -  9w(A +  C1S2  +  Di/32  +  Ei{%y-)2~ 

Flip i  ^  +  GrfW  ^)  +  2soV’(1)  ■ 

15t  =  +  #2^r  +  C2S2  +  D2P2  +  E2{^-)2- 

i^(1)  ^  +  G2^dA§~)  -  2s0^^, 


where  the  coefficients  are  lengthy  and  provided  in  Appendix  A.  This  system  of  equa¬ 
tions  is  linear  in  ( s2,/32 )  but  driven  by  nonlinear  functions  of  Vx  and  . 

The  steady  solution,  with  A  and  T  defined  in  the  appendix,  is 


(37) 


#2  (y)  =  A'i( 


cosh(\/  ErAy ) 
cosh(-\/  ErA ) 


i)  +  k2( 


cosh(\/  ErTy ) 
cosh(VrErr) 


1)  +  R\Er{y2  —  1), 


S2(y)  =  K3( 


cosh(-\/  ErAy ) 
cosh(VrErA) 


i)  +  -M 


cosh(\/  .ErTy) 
cosh(v^r) 


1)  +  S\Er{y2  —  1), 


where  the  coefficients  are  given  in  Appendix  B. 

We  denote  A2  =  l/{Erh)  and  p2  =  1  /(ErT);  A2  is  always  positive  and  concave  up 
as  a  function  of  9 ,  whereas  p2  is  monotonically  decreasing  and  may  change  sign  as  9 
varies.  For  example,  /12  goes  through  zero  at  a  critical  degree  of  anisotropy  9C  =  2.93 
for  N  =  6.  For  9  >  9C,  the  steady  state  becomes  highly  oscillatory;  we  show  in  the 
study  of  transient  solutions  that  this  behavior  coincides  with  the  onset  of  ill-posedness 
in  the  governing  system  of  equations.  The  behavior  of  9C  versus  concentration  N  can 
be  gleaned  from  Figure  3. 

Since  the  dominating  terms  in  the  order  parameters  near  y  =  ±1  are 


(38) 


evrErAy  eVErTy 

e\/~ErA  ’  gV^T  ’ 


the  order  parameters  have  a  boundary  layer  near  the  wall,  whose  width  is  proportional 
to 


(39) 


1  1 

'/ErA  ’  VEfT  ’ 


respectively.  These  are  the  penetration  depths  of  the  wall  layer  for  tangential  anchor¬ 
ing,  which  agrees  with  the  asymptotic  analysis  of  the  DMG  model  [14]  in  the  single 
elastic  constant  limit.  One  finds  the  order  parameters  are  coupled  with  anisotropic 
elasticity,  i.e. ,  the  orientational  distribution  is  strongly  biaxial  (birefringent  in  any 
plane).  For  both  rods  and  discs,  by  comparing  the  two  exponential  terms  in  and 

s2,  we  notice  that  the  boundary  layer  in  s2  is  governed  by  whereas  in  by 

so  that  their  scaling  behavior  is  incommensurate  with  the  leading  order  wall 
layer  scaling. 
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Table  1 

Steady  state  features  of  the  order  parameter  morphology  for  Couette  flow  with  tangential  an¬ 
choring  (BL  denotes  boundary  layer). 


FA/rods 

FA  /  discs 

T/rods 

T/discs 

S  -  so 

Concave  down 

Concave  up 

Concave  up 

Concave  down 

Sc  Concave  up  in  BL 

a 

Concave  up 

Concave  up 

Sc  Concave  down  in  BL 

Concave  down 

Concave  down 

-e 


x  10 


1 


(a) 


(b) 


xIO 


.4 

xIO 


W 

I 

w 


(e) 


Fig.  4.  The  steady-state  asymptotic  solutions  (3  (left  column)  and  s  —  so  (right  column)  in 
the  tumbling  regime  as  functions  of  (0,  y)  with  tangentially  anchored  boundary  conditions  in  weak 
Couette  flows,  (a)  and  (b)  depict  the  solution  for  rods  with  parameter  values  a  =  0.8 ,N  =  6 ,  De  = 
0.01,  Er  =  100.  (c)  and  (d)  depict  the  solution  for  discs  with  parameter  values  a  =  —0.8,  A  = 
6,  De  =  0.01,  Er  =  100. 


Table  1  tabulates  the  general  behavior  of  the  steady  states  in  the  regime  of  flow¬ 
aligning  and  tumbling  for  both  rods  and  discs,  in  which  FA  stands  for  flow-aligning, 
T  denotes  tumbling,  and  BL  denotes  a  boundary  layer.  The  concavity  switching 
phenomenon  between  tumbling  and  flow-aligning  materials,  noted  in  [14],  is  observed 
in  both  rods  and  discs.  The  steady  states  are  in  general  insensitive  to  changes  in  6.  In 
addition,  the  order  parameter  correction  S2  goes  through  a  similar  transition  for  rods 
versus  discs,  whereas  j3  has  the  same  concavity;  this  is  another  indication  of  strong 
biaxiality  for  anisotropic  elasticity.  Figure  4  depicts  a  typical  steady  state  asymptotic 
solution  of  (s2,/?2)  for  tumbling  rods  (a  =  0.8)  and  tumbling  discs  (a  =  —0.8). 
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2.2.4.  Transient  behavior  of  (S2,  (3?)  near  steady  states.  The  transient 
behavior  of  the  order  parameters  S2  and  /?2  obey: 


(40) 


dS2 

- 1  | 

(V2f 

dt  = 

9  Er  ' 

^  dy 2 

dy2 

d/?2  _ 

-1 

Udlsi 

+  B2d2J22 

dt 

9  Er 

l  dy2 

dy2 

+  C1S2  +  D1P2 


+  C2S2  +  D2P2 


with  zero  boundary  conditions. 

After  tedious  but  straightforward  Fourier  analysis,  the  growth  (er'  >  0)  or  decay 
(a'  <  0)  of  solutions  S2(y,t)  and  P2 (y, t),  due  to  the  time-dependent  factor  eat,  is 
determined  by  sgn(cr') 

<4  =  -5 [Ci  +  D2  -  Aik2  -  B2k2± 


V(Ci  +D2-  Aik2  -  B2k2)2  -  4 ((A2k2  -  C2){D1  -  B±k2)  +  (Ci  -  A1k2){D2  -  B2k2))}. 

(41) 

We  now  analyze  sgn(a'±)  to  deduce  stability.  For  long  waves  (|fc|  <<  1),  asymp¬ 
totic  formulae  can  be  derived: 

(42)  4  ~ -\[Ci  +  £>2  ±  ^{Ci  +  D2)2  +  4(C2Di  -  C1D2)]. 

Using  the  formulae  in  Appendix  B,  Ci,D2  >  0,(72  =  0,  which  implies  a'±  <  0.  On 
the  other  hand,  the  growth  rate  for  short  waves  (\k\  >>  1)  is  dominated  by 

(43)  cr±  ~  -[Ai  +  f?2  i  \/ (Ai  +  B2)2  +  4(A2l?i  —  AiB2)\k2. 


We  find  a±  <  0  only  for  sufficiently  small  |0|.  Otherwise,  A1B2  —  A2B1  may  be 
negative,  leading  to  a  positive  growth  rate  proportional  to  k2,  an  ill-posed  behavior. 
The  transition  to  ill-posedness,  A1B2  —  A2B1  =  0,  simplifies  dramatically  to  the 
condition  9  =  5so_2 ,  where  so  is  given  in  (8).  This  neutral  stability  curve  is  plotted 
in  Figure  3. 

Proposition  2.  The  steady  state  solution  (s2,/?2)  of  (36)  is  catastrophically 
unstable  if  and  only  if  the  degree  of  anisotropic  elasticity  satisfies 


(44) 


9  > 


6 

5s0  -  2' 


2.2.5.  Rheological  features  of  steady  structures.  The  shear  viscosity  (shear 
stress  divided  by  local  shear  rate)  at  the  plates  is  identical  to  the  averaged  shear 
viscosity  over  the  shear  cell;  it  is  a  nonzero  constant  at  O(De),  given  by 


(45) 


Vwall 


aa2(m^)[e(l-So)+3(W^-)] 
36(3+0(1  — s0)) 


+  |(Miso  +  3ry). 


It  can  be  readily  shown  that  r]waii  is  a  slowly  varying,  decreasing  function  of  the 
degree  of  anisotropic  elasticity  9  for  tumbling  rods  and  all  platelet  nematic  liquids; 
however,  rjwau  increases  versus  9  for  flow-aligning  rods.  Nonzero  first  normal  stress 
difference  Ni  and  second  normal  stress  difference  N?  show  up  at  order  0(De2)  and 
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are  given  in  Appendix  C.  These  non-Newtonian  effects  are  measurable  physically, 
e.g.,  Ni  >  0  corresponds  to  pushing  the  parallel  plates  apart  and  <  0  corresponds 
to  pulling  the  plates  together.  N\  is  a  linear  combination  of  a  and  fx 2  and  N2  is  a 
linear  combination  of  a  and  fx  1  + 112-  At  walls,  the  terms  containing  /j-1,2  drop  out  so 
that  N\  and  IV2  are  proportional  to  a: 

N1  =  [^(K3X2  +  K4fx2  +  2S) 

.  ( — 4a0+3so-l-16(i0SQ  —  12a+12asQ  —  12cl0sq )  /  -rr  \  2  1  ts~  ,,2  1  o  td\ 

+ - 108 E? -  (AlA 

+  427Er  (~4a0so2  —  4a0  +  24a0so3  —  6aso  +  9sq  +  18aso2  —  12a  —  16a0so)]a, 


No  =  \- 


bAEr 


(K3X2  +  /A4/x2  +  2  S) 


+  {&ae+12a-Z+&aso^l6as20-V2s0)s0  (^2  +  K^2  +  2  R) 

+  iifer(86's 02  +  26»  -  24(9s03  +  12s0  -  9s0  -  18as02  +  6  +  146>s0)]a, 

(46) 

where  G,  K ,  A,  S',  A,  p  and  Ki(i  =  1,  2,  3, 4)  are  given  in  Appendix  C. 

Figure  5  depicts  N\  and  N2  for  tumbling  rods  as  well  as  discs  as  functions  of  y  at 
some  parameter  values;  Figure  6  shows  N3  and  N2  for  flow-aligning  rods  and  discs. 
Table  3  lists  the  averaged  normal  stress  differences  calculated  in  four  representative 
cases.  We  summarize  the  noticeable  features  in  the  stress  differences  below. 

•  For  flow-aligning  rods,  N\  <  0  and  N2  >  0  across  the  gap;  their  signs  change 
in  plate  boundary  layers  for  small  degree  6  of  elastic  anisotropy. 

•  These  properties  reverse  for  flow-aligning  or  tumbling  discs,  which  may  ex¬ 
perience  sign  changes  in  N\  and  N2  in  the  middle  of  the  plate  gap  for  small 

\e\. 

•  For  tumbling  rods  and  all  discotic  NLCPs,  iVi  >  0  and  N2  <  0,  with  sign 
changes  for  rods  in  the  wall  boundary  layer  for  large  6 ,  and  sign  changes  for 
platelets  in  the  midgap  at  small  \6\. 

•  The  averages  across  the  gap  yield  Ni  >  0  and  N2  <  0  for  flow-aligning  rods, 
and  Ni  <  0  and  N2  >  0  in  all  other  cases. 

2.3.  Homeotropic  anchoring  (yto  —  f  )•  The  boundary  anchoring  condition 
affects  only  the  coefficients  of  the  governing  system  of  partial  differential  equations  at 
each  order.  The  structure  of  the  equations  at  O(De)  for  are  identical  to 

(18)  with  the  following  new  coefficients: 


(47) 


^  ~~  gw(s°  +  2)  (3  +  0(2sq  +  1)),  B  —  —  i(l  +  A l), 


C  —  isEr  [^(2so  +  1)Al  +  3(A L  +  1)], 


D  =  |(pi  So  +  377). 


The  nonzero  leading  order  solution  for  the  order  parameters,  primary  velocity 
component,  and  major  director  angle  are 


(48)  so  =  so,  Po  =  0,  v^\y)  =  y,  4>{1) (y)  =  MEr{y2  -  1), 

(49)  M  =  2  (s0  +  2)(3  +  9(2s0  +  1)) (1  +  Ai)' 
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Fig.  5.  The  normal  stress  differences  Ni  and  N2  in  tumbling  regimes  as  functions  of  y  at 
selected  values  of  6  with  tangentially  anchored  boundary  conditions  in  weak  shear,  (a)  Rods:  a  = 
0.8,  N  =  6 ,De  =  0.01,  Er  =  100,  a  =  10,  p.1  =  0.01, 77  =  0.02.  (b)  Discs:  a  =  -0.8,  N  =  6 ,  De  = 
0.01,  Er  =  100,  a  =  10,  m  =  -0.01,  y  =  0.02. 


It  is  easy  to  see  that  \M\  decreases  with  respect  to  9;  M  is  negative  for  flow-aligning 
discs  (o  <  0,  Al  <  —  1)  and  positive  for  all  other  cases.  Thus,  directors  wind  coun¬ 
terclockwise  from  the  lower  shearing  plate  to  the  midplane,  and  unwind  from  the 
midplane  to  the  upper  plate  for  flow-aligning  discs.  The  winding  reverses  direction  in 
the  other  cases. 

Analogous  to  tangential  anchoring,  steady  states  may  be  catastrophically  unstable 
for  discs  in  the  flow-aligning  regime  a  <  0,Xl  <  — 1  if  AD  —  BC  <  0,  yet  are  stable 
in  the  other  cases,  where 

AD  —  BC  =  2fEf  [(s  +  ty{ffiso  +  3?7)(0(2so  +  1)  +  3) 

(50) 


_l_aso  _|_  Aj,)(a0(2so  +  l)(so  +  2)  +  9sq(1  +  Al))]. 


Proposition  3.  The  steady-state  solution  is  stable  so  long  as 

(51)  9  >  9  = _ 12(s0  +  2)(/xis0  +  3??)  +  9aso(l  +  A l)2 _ 

4(so  +  2)(2so  +  l)(/iiSo  +  377)  +  aaso(l  +  Al)(1  +  2so)(2  +  So) 

The  governing  system  of  equations  for  the  order  parameters  (s,  (3)  at  order  0(De2) 
is  of  the  same  form  as  tangential  anchoring,  with  different  coefficients  given  in  Ap¬ 
pendix  C. 
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Fig.  6.  The  normal  stress  differences  Ni  and  N2  in  flowing-aligning  regimes  as  functions  of 
y  at  selected  values  of  6  with  tangentially  anchored  boundary  conditions  in  weak  shear,  (a)  Rods: 
a  =  0.9,  N  =  6,  De  =  0.01,  Er  =  100,  a  =  10,  /11  =  0.01,  rj  =  0.02.  (b)  Discs:  a  =  -0.9,  N  =  6,  De  = 
0.01,  Er  =  100,  a  =  10,  p.1  =  -0.01,  rj  =  0.02. 


The  steady  solutions  are 


(52) 


fc(y)  =  Kn( 


cosh(\/  ErAy ) 
cos  h(VErA) 


l)  +  R2Er(y2 -1), 


S2(y)  =  K21( 


cosh(\/  ErAy ) 
cosh(\/  ErA ) 


1)  +  A'22( 


cosh(\/  ErTy ) 
cosh(\/  ErT) 


1)  +  S2Er(y2  -  1). 


In  this  solution,  y2  (defined  earlier)  changes  sign  only  as  9  varies  below  a  threshold 
value  0d  for  platelets.  Again ,  the  change  of  sign  in  r2  coincides  with  ill-posedness  in 
the  governing  system  of  equations. 

Notice  that  in  s2,  there  are  two  cosh  terms,  whereas  there  is  only  one  in  f}2 .  For 
rods,  the  second  term  in  s2  dominates  in  the  boundary  layer  while  the  first  term 
dominates  for  discs.  Table  2  tabulates  features  of  the  steady  states.  Compared  with 
results  above  for  tangential  anchoring,  the  steady  states  with  normal  anchoring  are 
more  sensitive  to  the  degree  9  of  elastic  anisotropy.  The  order  parameter  variation 
versus  9  decreases  for  rods  and  increases  for  platelets.  The  solution  profiles  switch 
concavity  in  the  boundary  layer  for  flow-aligning  versus  tumbling  discotics,  whereas 
the  concavity  remains  the  same  for  rods.  Figure  7  depicts  typical  steady  solutions  for 
tumbling  rods  and  discs. 
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Table  2 

Steady-state  features  of  the  order  parameter  morphology  for  Couette  flow  with  normal  anchoring. 


FA/rods 

FA/discs 

T  /  rods 

T  /  discs 

s  —  so 

Concave  down 

Concave  down 
&  Concave  up  in  BL 

Concave  down 
&  Concave  up  in  BL 

Concave  up 

&  Concave  down  in  BL 

0 

Concave  up 

Concave  up 

Concave  up 

&  Concave  down  in  BL 

Concave  down 

Fig.  7.  The  steady-state  asymptotic  solution  as  functions  of  ( 0,y )  in  tumbling  regime  with 
normally  anchored  boundary  condition  in  weak  shear,  (a)  and  (b)  depict  the  solution  for  rods  with 
parameter  values  a  =  0.8,  N  =  6,  De  =  0.01,  Er  =  100.  (c)  and  (d)  depict  the  solution  for  discs  with 
parameter  values  a  =  —0.8,  N  =  6,  De  =  0.01,  Er  =  100. 


Again  the  governing  system  can  become  ill-posed.  The  growth  rate  formulae  for 
the  steady  states  are  identical  to  (41),  but  with  new  coefficients.  Notice  that  the 
discriminant  A\B2  —  A2B i  =  4(1  —  so)2(l  +  2so)[3  +  0(1  +  2s0)][3  +  0(1  +  4so)].  For 
rods,  6  >  0  and  AiB2  —  A2B\  >  0  indicating  stability.  For  discs  and  0  <  —  1+4Sq  ,  the 
discriminant  is  negative  so  that  the  steady  state  is  unstable. 

Proposition  4.  The  system  for  the  two  order  parameters  (s2,/32)  is  locally  ill- 
posed  if  and  only  if  0  <  Og,  where 


(53) 


0d 


3 

1  +  4s0 ' 


The  shear  viscosity  at  the  walls  in  this  case  is  given  by 


(54) 


TJwall  — 


dy 


«a2(1+Ay)(2e(2so  +  1)+3(1+5y)) 
36(6»(2s0+l)+3) 


+  \(ld\SQ  +  377). 
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Fig.  8.  The  normal  stress  differences  Ni  and  N?  as  functions  of  y  in  tumbling  regime  at 
selected  values  of  0  with  normally  anchored  boundary  conditions  in  weak  shear.  The  parameter 
values  are  a  =  ±0.8,  N  =  6,  De  =  0.01,  Er  =  100,  a  =  10,  pi  =  ±0.01,  p2  =  0.02  for  rods  and  discs, 
respectively. 


Similar  to  the  case  of  tangential  anchoring,  the  viscosity  decreases  with  respect  to  9 
for  rods  and  flow-aligning  discs,  but  increases  with  respect  to  9  for  tumbling  discs. 
The  leading  order  normal  stress  differences  again  show  up  at  0(De 2)  and  are  given 
in  Appendix  D. 

Figures  8  and  9  depict  representative  plots  of  the  first  and  second  normal  stress 
differences  in  tumbling  and  flow-aligning  regimes,  respectively,  and  demonstrate  the 
following  properties: 

•  Rods  are  more  sensitive  to  the  variation  in  degree  9  of  elastic  anisotropy  than 
discs. 

•  For  tumbling  rods,  N\  >  0  and  N2  <  0  except  there  may  be  a  sign  change  in 
a  boundary  layer  near  the  wall  at  large  9.  For  tumbling  discs,  N\  and  IV2  are 
both  negative  except  that  is  positive  in  a  boundary  layer  near  the  wall. 

•  The  behavior  of  the  normal  stress  differences  does  not  change  much  for  flow¬ 
aligning  rods  compared  with  tumbling  rods.  For  flow-aligning  discs,  N\  >  0 
and  IV2  >  0  except  that  N2  is  negative  in  boundary  layer  near  the  wall. 

•  The  average  values  across  the  gap  obey  N\  >  0  and  IV2  <  0  for  rods,  but  they 
may  change  signs  for  discs. 

The  averaged  normal  stress  differences  are  tabulated  in  Table  3. 

In  summary,  the  salient  predictions  from  this  analysis  are: 
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Fig.  9.  The  normal  stress  differences  Ah  and  Ah  as  functions  of  y  in  flowing- aligning  regime 
at  selected  values  of  6  with  normally  anchored  boundary  conditions  in  weak  shear.  The  parameter 
values  are  a  =  ±0.9,  N  =  6,  De  =  0.01,  Er  =  100,  a  =  10,  pi  =  ±0.01,  p2  =  0.02  for  rods  and  discs, 
respectively. 


Table  3 

The  averaged  normal  stress  differences  (Couet.te). 


FA/rods 

FA/discs 

T  /  rods 

T  /  discs 

Tangential 

Ah  <  0,  N2  >  0 

N\  >  0,  N2  <  0 

Ah  >  0,  Ah  <  0 

N\  >  0,  N2  <  0 

Normal 

Ah  >  0,  Ah  <  0 

Ah  >  0 

Ah  >  0,  Ah  <  0 

iVi  <  0,  N2  <  0 

•  The  major  director  winds  counterclockwise  from  the  bottom  to  top  plates, 
for  both  flow-aligning  rods  in  tangential  anchoring  and  flow-aligning  discs  in 
homeotropic  anchoring.  Remarkably,  the  principal  orientation  axis  rotates 
clockwise  if  the  nematic  polymer  tumbles  in  weak  shear  rather  than  flow 
aligns.  The  magnitude  of  winding  of  the  orientation  axis,  which  sets  the 
number  of  bands  of  nematic  distortion,  reduces  with  the  degree  of  elastic 
anisotropy  6. 

•  The  order  parameters  are  relatively  insensitive  to  the  degree  of  elastic  anis¬ 
otropy  in  tangential  anchoring,  and  more  sensitive  in  normal  anchoring. 

•  Ill-posedness  may  occur  within  each  order  of  asymptotic  equations  depending 
on  the  values  of  the  parameters,  although  the  full  equations  are  well  posed. 
This  transition  implies  a  breakdown  in  the  asymptotic  ordering  which  allows 
explicit  solution  and  scaling  properties,  and  suggests  a  physical  transition 
away  from  these  asymptotic  structures. 
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v=0,  Q=Q0 


Fig.  10.  The  geometry  of  the  plane  Poiseuille  flow.  A  pressure  gradient  ^  =  —  De2  is 
imposed  across  the  channel.  At  the  bounding  surfaces,  the  orientation  tensor  is  assumed  to  equal  to 
its  equilibrium  value. 


•  The  averaged  shear  viscosity  varies  weakly  with  respect  to  the  strength  of 
anisotropic  elasticity. 

•  The  averaged  normal  stress  differences  may  take  on  all  possible  signs  depend¬ 
ing  on  the  parameter  regime. 

3.  Spatial  structures  in  weak  Poiseuille  flows.  In  this  section,  we  study 
steady  structures  in  the  direction  of  the  velocity  gradient  under  an  imposed,  small 
pressure  gradient  ^  in  plane  Poiseuille  flow.  Figure  10  depicts  the  cross  section  of  the 
plane  Poiseuille  flow  on  the  (x,  y)  plane.  The  boundary  condition  for  the  orientation 
tensor  is  identical  to  that  used  for  shear  flows  while  the  velocity  boundary  condition  is 
no-slip  v(±l)  =  0.  As  before,  we  use  h  and  to  =  tn  =  jj-  as  the  characteristic  length 
and  time  scale,  respectively.  We  adopt  the  same  dimensionless  symbols  used  in  weak 
plane  Couette  and  assume  =  —e  =  — De 2  in  the  dimensionless  form,  where  the 
Deborah  number  is  defined  by 


(55) 


De  = 


e  =  De2. 


We  seek  asymptotic  solutions  in  powers  of  e.  The  momentum  equation  yields  at  order 
0(e): 


(56) 


a, ,(!)  Oi-U) 

OVx  _  _ 1  U  xy 

dt  1  dy  • 


The  other  governing  equations  are  identical  to  those  derived  for  plane  Couette  flows. 
Hence,  the  transient  solution  and  the  stability  of  steady  states  are  identical  to  the 
corresponding  problems  in  plane  Couette  flows.  We  will  not  repeat  them  here;  in¬ 
stead,  we  only  present  the  asymptotic  steady  states  with  respect  to  the  two  anchoring 
conditions. 

3.1.  Tangential  anchoring  (tpo  —  0).  The  steady  solutions  up  to  order  0(e) 
are 


(57)  s0  =  s0,  /3o  =  0,  u£1}(y)  =  Hx(l  -  y2),  ^(1)(f/)  =  H2Ery(l  -  y2), 
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where  H\  and  H2  are  given  in  Appendix  E.  The  positivity  of  Hi  coincides  with  the 
stability  of  the  steady  state,  giving  rise  to  a  parabolic  velocity  profile.  H2  is  positive 
for  flow  aligning  rods  and  negative  otherwise  in  stable  steady  states.  H2  behaves  more 
or  less  like  the  diagnostic  M  in  the  plane  Couette  flow.  Notice  that  ip^  is  an  odd 
function  of  y  leading  to  an  asymmetric  major  director  pattern,  known  as  a  chevron 
pattern,  with  respect  to  the  midplane  [6,  2,  3]. 

For  flow-aligning  rods  in  stable  steady  states,  Hi  ( H2 )  decreases  (increases)  with 
respect  to  9,  and  increases  (decreases)  with  respect  to  9  in  all  other  cases.  The 
rotational  pattern  of  the  major  director  (a  function  of  ip^)  is  dictated  by  A l-  For 
flow-aligning  rods,  Xl  >  1,  the  major  director  rotates  counterclockwise  from  the  lower 
plate  to  the  of  the  shear  cell  and  then  reverses  its  rotation  to  the  midplane.  The 
orientation  pattern  in  the  top  half  of  the  cell  is  the  mirror  image  of  that  in  the  lower 
half.  The  rotation  reverses  for  the  other  cases  where  Xl  <  1. 

The  steady  solutions  of  the  order  parameters  at  order  0(e2)  are  given  by 


#2  (y)  =  Kn( 


cosh(\/  ErAy ) 
cosh(\/  ErA ) 


1  )  +  k12{ 


cosh(\/ ErYy) 
cosh(-\/  ErT) 


1)  +  Rn(y4-1)  +  Ri2(y2 -1), 


S2(y)  =  k21( 
(58) 


cosh(\/  ErAy ) 
cosh(\/  ErA ) 


1)  +  k22( 


cosh(\/  ErY  y) 
cosh(-\/  ErY) 


l)  +  R2i(y4-l)  +  R22(y2-l), 


where  the  coefficients  are  given  in  Appendix  E.  The  order  parameters  behave  like  a 
quartic  polynomial  with  respect  to  y  in  most  part  of  the  cell  except  at  the  boundary 
layers  near  the  plates.  The  velocity  v^\  angle  variable  ip ^  and  the  order  parameter 
/?2  are  insensitive  to  the  variation  of  9.  The  sensitivity  of  the  order  parameter  cor¬ 
rection  s2  is  the  most  pronounced  at  9  =  0,  i.e.,  in  the  one-constant  approximation. 
Figure  11  depicts  typical  steady  solutions  for  tumbling  rods.  Table  4  tabulates  all  the 
steady  state  behavior: 

•  The  thickness  of  the  boundary  layers  in  this  flow  are  narrower  than  those 
in  the  weak  plane  Couette  flows,  suggesting  a  mollifying  effect  of  stronger 
velocity  gradients  near  the  walls. 

•  The  two  order  parameter  corrections  (at  small  |d|)  and  the  angle  variable 
change  their  signs  and  concavity  in  the  tumbling  versus  flow-aligning  regime, 
but  only  for  rods;  this  predicts  a  profile  concavity  flip  in  the  focusing  and 
defocusing  of  the  orientation  distribution  occurs  as  rods  pass  through  the 
flow-aligning  to  tumbling  transition;  the  profile  of  (32  is  either  W-shaped  or 
M-shaped. 

•  The  velocity  profile  is  concave  down. 

•  The  parameter  s  —  so  is  very  sensitive  around  9  =  0  for  both  rods  and  discs. 

•  The  angle  profile  is  a  rotated-S  shape. 

Figure  11  depicts  typical  steady  solutions  for  tumbling  discs  as  functions  of  {9,y). 

The  shear  viscosity  is  given  by 


(59)  llapp  —  ~^TT  —  —  TShT  K«o  +  2)(0(1  —  so)  +  3(1  —  y^))]  +  |(miSo  +  3r/), 

dy 


a  decreasing  function  with  respect  to  9  for  discs  and  tumbling  rods  but  an  increasing 
function  for  flow-aligning  rods. 

The  normal  stress  differences  in  this  case  are  given  in  Appendix  E.  Unlike 
weak  Couette  flows,  they  are  rational  functions  of  yi,^2  and  a.  Figures  12  and  13 
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(C) 


Fig.  11.  Steady  solutions  of  /3,  s— so, v  andif—tfo  as  functions  of  (6 ,  y)  in  the  regime  of  tumbling 
rods  with  tangentially  anchored  boundary  conditions  in  plane  Poiseuille  flow.  The  parameter  values 
are  N  =  6,  a  =  0.8,  De  =  0.01,  Er  =  100,  pi  =  0.1, 77  =  0.2,  a  =  10. 


depict  some  representative  plots  of  the  normal  stress  differences  for  tumbling  and 
flow-aligning  nematics,  respectively.  In  summary,  they  have  the  following  properties: 

•  For  flow-aligning  rods,  Ni  is  positive  and  IV2  is  negative.  The  signs  are 
reversed  for  flow-aligning  discs. 

•  For  tumbling  rods  and  discs,  Ni  is  negative,  but  N2  is  positive. 

•  In  both  stable  tumbling  and  flow-aligning  regimes,  for  rods  and  discs,  the  ab¬ 
solute  values  of  iVi  and  N2  increase  and  decrease,  respectively,  as  a  increases. 
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Table  4 

Steady  states  (Poiseuille)  in  tangential  anchoring. 


FA/rods 

FA/discs 

T/rods 

T  /  discs 

s  -  So 

Concave  up 

Concave  up 

Concave  down  at  small  6 
&  Concave  up  at  large  6 

Concave  up 

p 

W-shape 

M-shape 

M-shape 

M-shape 

-g- 

1 

■g- 

0 

s 

Rotated-S 

Rotated-S 

Rotated-S 

Fig.  12.  The  normal  stress  differences  N\  and  N2  as  functions  of  y  in  tumbling  regime  at 
selected  values  of  6  with  tangentially  anchored  boundary  conditions  in  plane  Poiseuille  flows.  The 
other  parameter  values  are  N  =  6,De  =  0.01 ,  Er  =  100,  a  =  10,^2  =  0.02,  for  rods  a  =  0.8,  pi  = 
0.01,  and  for  discs,  a  =  — 0.8,p.i  =  —0.01. 


•  The  gap  averages  satisfy  N±  >  0  and  N2  <  0  for  flow-aligning  rods  while 
N\  <  0,  N2  >  0  in  all  other  regimes. 

The  behavior  of  the  averaged  normal  stress  differences  is  tabulated  in  Table  6. 


3.2.  Homeotropic  anchoring  (xp  —  ^ ) .  The  steady  solutions  up  to  order  O(e) 
are  given  by  (57)  with  new  H\  and  H2  given  in  Appendix  F.  As  in  the  tangential 
anchoring  case,  Hi  is  positive  in  all  stable  steady  states.  H2  is  negative  for  flow- 
aligning  discs  and  positive  otherwise  for  other  stable  steady  states.  For  flow-aligning 
discs  in  stable  steady  states,  Hi  (H2)  decreases  (increases)  with  respect  to  9  for  flow¬ 
aligning  discs,  yet  increases  (decreases)  with  respect  to  9  for  all  other  stable  steady 
states. 
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Fig.  13.  The  normal  stress  differences  Ni  and  N2  as  functions  of  y  in  flow- aligning  regime  at 
some  values  of  9  with  tangentially  anchored  boundary  conditions  in  plane  Poiseuille  flows.  The  other 
parameter  values  are  N  =  6,  De  =  0.01 ,Er  =  100,  a  =  10,  p2  =  0.02,  for  rods  a  =  0.9,  pi  =  0.01, 
and  for  discs,  a  =  —  0.9,  pi  =  —0.01. 


Table  5 

Steady  states  (Poiseuille)  in  homeotropic  anchoring. 


FA/rods 

FA  /  discs 

T / rods 

T  /  discs 

S  -  so 

W-shape 

M-shape 

W-shape 

W-shape 

(3 

W-shape 

M-shape 

W-shape 

W-shape 

-g- 

1 

O 

Rotated-S 

s 

Rotated-S 

Rotated-S 

The  steady  solutions  of  the  order  parameters  at  0(e2)  are 


#2(2/)  =  A"n( 


cosh(\/  ErAy) 

cos  h(^ErA) 


1)  +  R1(y^-1)  +  S1(y2-1), 


(60) 


S2(y)  =  K 21  ( 


cosh(\/  ErAy ) 
cosh(-\/  ErA) 


1) 


-i)  +  - ')  +  W  -  d- 

We  summarize  the  features  of  stable  steady  states  in  Table  5: 

•  The  two  order  parameters  and  the  angle  parameter  change  their  signs  and 
shapes  in  tumbling  and  flow-aligning  regime  for  discs  but  not  for  rods,  indi¬ 
cating  the  solutions  are  more  sensitive  for  platelet  molecules  than  for  rods. 
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Fig.  14.  Steady  solutions  /3,  s  —  .so.  v  and  %b  —  V'Q  for  tumbling  rods  as  functions  of  ( 9,y ) 
with  normally  anchored  boundary  conditions  in  plane  Poiseuille  flows.  The  parameter  values  are 
a  =  0.8,  N  =  6,  De  =  0.01,  Er  =  100,  mi  =  0.01,  r?  =  0.02,  a  =  1. 


•  As  the  anisotropic  elasticity  enhances,  the  order  parameter  variations  and  the 
angle  variation  reduces  for  rods  yet  amplifies  for  discs. 

•  The  velocity  profile  has  fixed  concavity  in  all  regimes.  As  the  anisotropic 
elasticity  increases,  the  velocity  increases  for  rods. 

•  The  orientational  variables  decrease  with  respect  to  0  while  the  flow  variable 
v  increases. 

•  Again,  the  thickness  of  the  boundary  layers  are  smaller  compared  to  weak 
plane  Couette  flows. 

Figure  14  depicts  a  typical  steady  solution  for  tumbling  discs  as  functions  of  {0,y). 
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Fig.  15.  The  normal  stress  differences  N\  and  N2  as  functions  of  y  in  tumbling  regime  at 
selected  values  of  6  with  normal  anchoring  boundary  conditions  in  plane  Poiseuille  flows.  The 
parameter  values  are  N  =  6,  De  =  0.01,  Er  =  100,  a  =  10,,  a  =  0.8,  pi  =  0.01,  ^2  =  0.02  for  rods, 
and  a  =  — 0.8,/ii  =  —  0.01,  P2  —  0.02  for  discs. 


The  wall  shear  viscosity  is  given  by 

(61)  Vwall  =  =  18J5VHV  (a(2  +  so)(0(2so  +  1)  +  3(1  +  yjr)))  +  |(miso  +  3??), 

dy 

which  decays  with  respect  to  9  for  all  rods  and  flow-aligning  discs,  but  increases  for 
tumbling  discs.  The  first  and  second  normal  stress  differences  in  this  case  are  given  in 
Appendix  F.  Figures  15  and  16  depict  the  normal  stress  differences  for  flow-aligning 
and  tumbling  nematics,  respectively.  In  summary,  they  exhibit  the  following  features: 

•  For  flow-aligning  rods  and  disks,  N\  and  N2  are  negative  except  for  a  small 
region  at  the  midplane,  where  some  of  the  first  normal  stress  difference  may 
be  positive. 

•  For  tumbling  rods,  N\  and  IV2  are  negative.  For  tumbling  discs  (— 1  <  < 

0),  Ni  and  N2  are  negative  except  for  a  tiny  region  at  the  midplane  at  small 
9.  For  tumbling  discs,  the  behavior  reverses  completely,  i.e.,  the  normal  stress 
differences  are  positive  except  for  a  small  region  at  the  midplane. 

•  In  both  stable  tumbling  and  flow-aligning  regimes,  for  rods  and  discs,  the 
absolute  values  of  Ni  and  N2  will  decrease  and  increase  as  a  increases,  re¬ 
spectively. 

•  The  gap  averages  obey  Ni  >  0  and  N2  >  0  for  tumbling  rods  and  both 
become  negative  in  all  other  regimes. 

The  behavior  of  averaged  normal  stress  differences  is  tabulated  in  Table  6. 
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Fig.  16.  The  normal  stress  differences  N\  and  N2  as  functions  of  y  in  flow- aligning  regime 
at  selected  values  of  6  with  normal  anchoring  boundary  conditions  in  plane  Poiseuille  flows.  The 
parameter  values  are  N  =  6,  De  =  0.01,  Er  =  100,  a  =  10,  a  =  0.9,  pi  =  0.01, 112  =  0.02  for  rods, 
and  N  =  6,  De  =  0.01,  Er  =  100,  a  =  1,  a  =  —0.9,  pi  =  —0.01,  P2  —  0.02  for  discs. 


Table  6 

The  normal  stress  differences  (Poiseuille). 


FA/rods 

FA/discs 

T  /  rods 

T  /  discs 

Tangential 

Ah  >  0,  Ah  >  0 

Ah  <  0,  Ah  >  0 

Ah  <  0,  Ah  >  0 

Ni  <  0,  N2  >  0 

Normal 

Ah  <  0,  Ah  <  0 

Ni  <  0,  N2  <  0 

Ah  <  0,  Ah  <  0 

Ni  >  0,  N2  >  0 

4.  Conclusion.  We  have  derived  explicit  asymptotic  structures  for  weakly 
sheared  nematic  polymers  in  both  plate-driven  and  pressure-driven  experimental  con¬ 
ditions.  The  goal  of  this  analysis  is  to  predict  scaling  properties  in  the  orienta¬ 
tional  distribution  of  the  rigid  rod  ensemble  from  the  strong  elasticity,  weak  flow 
regime,  which  then  guide  numerical  continuation  studies  of  heterogeneous  films  and 
molds  across  a  multi-parameter  space  of  material  properties  and  processing  conditions. 
We  have  explored  the  effect  of  anisotropic  elasticity  for  both  flow  conditions,  using 
a  second-moment  model  for  the  orientational  distribution  derived  from  Doi  Hess- 
Marrucci-Greco  kinetic  theory.  These  results  extend  previous  work  of  the  authors  for 
steady  structures  with  equal  elasticity  constants  in  plate-driven  flow  in  several  ways: 
anisotropic  elasticity,  pressure-driven  flows,  and  transient  asymptotic  equations  fol¬ 
lowed  by  stability  predictions  within  the  asymptotic  ordering  of  the  flow-nematic  sys¬ 
tem.  The  leading  order  flow  structure  is  simple  linear  shear  versus  Poiseuille  profiles 
for  the  respective  driving  conditions,  with  elastic  hydrodynamic  feedback  contribu¬ 
tions  characterized  at  next  order.  These  results  confirm  the  consistency  of  imposing 
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the  kinematics  (and  thereby  suppressing  flow  feedback)  in  special  asymptotic  parame¬ 
ter  regimes,  yet  also  predict  breakdown  of  this  decoupling  of  the  momentum  equation 
when  either  of  several  conditions  is  relaxed:  weak  flow,  strong  elasticity,  or  sufficiently 
isotropic  elasticity. 

The  orientational  structures  for  both  flow  conditions  convey  scaling  properties  of 
nematic  (director  dominated)  elastic  distortions  as  well  as  molecular  elasticity  (domi¬ 
nated  by  focusing  or  defocusing  of  the  orientational  distribution) .  The  structure  scal¬ 
ing  laws  are  similar  for  plane  Couette  and  Poiseuille  flows,  with  plate  boundary  layer 
thicknesses  proportional  to  1/ \J  Er  and  nonuniform  structures  spanning  the  plates 
with  mean  lengthscale  proportional  to  1/ Er.  The  prefactors  of  the  structures  capture 
the  roles  of  material  properties:  flow-aligning  versus  tumbling  nematics,  degree  of 
anisotropy  in  the  elasticity  potential,  strength  of  the  short-range  nematic  potential, 
and  molecular  aspect  ratio.  These  subtleties  are  detailed  in  the  body  of  the  paper, 
where  the  amplitude  of  structure  variations,  convexity  of  profiles,  and  stability  of  the 
steady  structures  all  depend  strongly  on  these  molecular  parameters  as  well  as  plate 
anchoring  conditions. 

The  particular  results  are  less  important  than  the  overall  insight  into  the  sen¬ 
sitivity,  flow-nematic  feedback  and  processing-generated  structures  on  material  and 
device  properties.  The  instability  within  the  asymptotic  equations  is  catastrophic, 
similar  to  backward  heat  flow  instabilities,  so  there  is  no  mistaking  the  breakdown  of 
these  steady  profiles  in  the  weak  flow  model  system.  The  role  of  anisotropic  elastic¬ 
ity  is  shown  to  be  greater  for  normal  versus  tangential  plate  anchoring,  and  greater 
in  pressure-driven  than  plate-driven  flows.  In  Poiseuille  flow,  the  transition  to  catas¬ 
trophic  instability  coincides  with  a  rapidly  growing  midplane  axial  velocity,  confirming 
a  breakdown  in  the  asymptotic  analysis.  Finally,  anisotropic  elasticity  is  shown  to 
contribute  to  either  shear  thinning  or  shear  thickening  behavior  as  other  parameters 
are  modified,  and  signs  of  normal  stress  differences  (which  determine  whether  the 
plates  are  pushed  or  pulled  by  the  stresses  generated  between)  are  likewise  sensitive 
to  various  material  parameters. 

Appendix  A.  The  coefficients  in  the  second-order  equations  in  tangen¬ 
tial  anchoring. 

Ai  =  (1  —  s0) [ — 6(1  +  2s0)  +  (8s§  -  so  -  2)9],  B1  =  (1  -  s0)s0(6  -  7s06), 

Ci  =  18ErNs0(-l  +  4s0),  D1  =  36ErNs0{l  -  s0), 

Ei  =  2s0(l  —  so)(l  +  3so)[6  +  (2  —  5so)0], 

(62)  Fi  =  2so[3(2  +  so)  +  2(3sq  +  so  —  1)0],  Gi  =  9£Vso( — 1  —  o  +  2aso), 

A2  =  — 2so0(l  -  s0),  B2  =  (1  -  s0)[— 6  +  (3s0  —  2)0],  C2  =  0, 

D2  =  54 ErNs0,  E2  =  2 s0(s0  -  1)[6  +  (2  -  50)], 

F2  =  2s0[3(2  -  s0)  +  (so  +  2)(s0  -  1)0],  C2  =  9Ers0(l  -  a). 

Appendix  B.  The  coefficients  in  the  second-order  equations  in  normal 
anchoring. 

A\  =  —2(1  —  so)(2so  +  1)[3  +  (4so  +  1)0], 

B\  =  2s0(l  —  s0)  [3  +  2(2s0  +  1)0], 

Ci  =  18£VlVso(— 1  +  4s0),  D1  =  36£WVs0(l  -  s0), 

,  ,  Ei  =  4s0(l  -  s0)[3(l  +  3s0)  -  (2s0  +  l)(6s0  +  1)0)], 

1  j  Fi  =  2s0[3(2  +  s0)  -  (6sg  -  5sg  -  8s0  -  2)0],  Gi  =  9£Vs0(-l  +  a  -  2 as0), 

A2  =0,  B2  =  2(s0  -  1)[3  +  (2s0  +  1)0],  C2  =  0,  D2  =  54 ErNs0, 

E2  =  4so[3(— 1  +  so)  +  (2s§  +  3so  —  1)0], 

F2  =  2s0[— 3(2  +  s0)  +  (2s0  +  l)(s0  -  2)0],  G2  =  9Crs0(l  +  a). 
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Appendix  C.  The  coefficients  of  the  steady  solutions  for  tangential 
anchoring  in  weak  plane  Couette  flows. 

A ,fj,  =  \/Er{18[Ns0(so  -  l)((16s02  -  5s0  -  2 )9  -  6(5s0  +  1))t 

^((eSso4  +  49s03  -  24s02  -  20so  +  8)92  -  6(8s03  +  21s02  +  6s0  -  8)6>+ 

18(s0  +  2)2)1/2So(i  _  s0)N]/[{so  -  l)2(2s0  +  l)(50*o  -2 9-  6 )(9s0  -29-  6)]}i, 

T-,  _  (aso— 3so+2a)(so_ 1)[(— 5so+2asQ+3aso+2so— 2a)0— 3(aso— 2so+2a)] 

Ul  ~  4iVs2(so+2)2(0so-0-3)2  ’ 

Si  =  (1  —  so)[(— 88sq  +  34aso  +  77asg  —  31sg  +  8aso  +  20so  —  20a)6>  —  6(7asg  —  14sg 
—  IOsq  +  19aso  +  10a)]  (aso  —  3sq  +  2a)/[8s§lV(— 1  +  4sq)(sq  +  2)2(9sq  —  9  —  3)2], 


A  =  ^,  T  =^,R=  RiEr,  S  =  SiEr,  K  =  , 


VEr 


4so(so+2)[(sq  —  1)6— 3]  ' 


A  =  -(s0  -  l)2(2s0  +  l)(56»s0  -29-  6 )(9s0  -29-  6), 

B  =  36ErNso(so  -  l)[(16s§  -  5s0  -  2)9  -  6(5 s0  +  1)],  C  =  972(£rAs0)2(l  -  4s0), 

rp  _  2A1(4K2E2+2K2F2+KG2)-2A2(4K2E1+2K2F1+KG1)-C1(2K2F2+KG2)-2BR 
1  ~  C  > 

Ki  =  i  [«24ifc«2Mi  +  fj2 (R  +  T)  -  2R] ,  K3  =  , 


K,  = 


A2  /t2  L  A  r  nyiv-rs-j  “3  -  42A2 

[4K2A1E2—4K2A24E1  +A2(j^  +  r)-2  R],  K4  =  -B2£+?2K2, 


-\2+fl2 1 


1U  _  Ga  ( K3\2  cosh  (A  y)  .  K4[j2  cosh  (/j,y)  ,  9  rA 

1  54 Er  y  cosh(A)  '  cosh  fi  '  ) 


aNso(l—4  so)a 
3 


K:> 


( 


cosh(A  y ) 
cosh(A) 


-i  + 


k4  ( 


cosh (/j,  y ) 
cosh (y) 


-i)  +  s  (y2~  i) 


■  (— 4a0+3so  +  16a0Sg  —  12a+12asQ  —  12a0so)a  (  K\X2  cosh(Ai/)  .  K2fi2  cosh (fjy)  ,  cy  td\ 
'  108  Er  y  cosh(A)  '  cosh  (fi)  '  J 


Afs0(l+2  s 0)  aOL  ^ 


cosh(A  y) 
cosh(A) 


-1  )+k2  (- 


:osh(/x  y) 
cosh  (//) 


-1)  +R{y2-  1) 


+  4^7J“a  [4a(so  —  l)(2so  +  l)(3so  +  1)9  —  3(2aso  —  3sg  —  6aso2  +  4 a]y2 

[(2s0  +  l)(3so2  -  5s0  -  4))9  -  6 (s0  +  2)](y2  -  1)  +  2/x2  Ks02  ( y 2  -  l), 

at  _  _  Ga  ( K3\2  cosh(Ay)  .  Kgj2  cosh ()iy)  ,  9  q\ 

2  54 Er  cosh(A)  '  cosh /r  '  J 


,  aJVso(-l+4s0)o: 
3 


K* 


sh(A  y) 
3sh(A) 


(8a0+12a— 3+8aso  — IhaSg  —  12so)sqQ: 


10  8Er 


2Nso(l—so)act 

3 


Ki 


cosh(A  y) 
cosh(A) 


-  l)  +I<4  (; 

(& 

)+M 


sh(/f  y) 


coshoir  -  !J  +  5  (y2  - !) 

^  .  2 _ T_/..  _.\ 

2  R 


A2  cosh(Ay)  .  K2IA  cosh (fiy) 
cosh(A)  '  cosh  (/i) 


cosh  (/x  y) 
cosh  (/x) 


-i)+i?(j/2-i) 


+  4 27 Er*  [~2a(so  -  l)(2so  +  l)(6s0  +  1)9  +  3(4as0  -  3s0  -  6as02  +  2 a)}y2 
+  2°27 Er  a  [—  (^s°  4”  l)(3so2  —  4so  —  2 ))9  +  3(so  +  2 )](y2  —  1) 

-2  (th  +  92)Ks02  ( y 2  -  1)  , 

G  =  2a9  +  4a0Sg  +  6aso  —  3so  —  12asg  —  16a0Sg  +  lOa0so  +  6a. 

(64) 
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Appendix  D.  The  coefficients  of  the  steady  solutions  in  normal  anchor¬ 
ing  in  weak  plane  Couette  flows. 


a  =  v^[(1_so)^:°o+3+e)]1/2,  n  =  VE-ri^ 


9Ns0(  —  1+4sq) 


(l-so)(2so  +  l)(40so+3+e)J 


11/2  *  _  _A_  p  _  _V_ 


^/Er, 


p  9(1  — so)(25o+qso+2q)(qso+38q+2q)  p  p  j?.  p-  3(qgo+3go+2a)  p 

n  ~  4N(s0+2)2(29s0+3+6)sl  rj' '  1X2  ~Er'  ~  4s0  (s0+2)  (20so+e+3)  ^ '  ’ 

T  =  3( — 1  +  So)[2(2so  +  l)(so  —  l)(aso  T  2so  +  2 a)0  +  ErNsoa^sg  +  2) 

+6( — 1  +  so)(2so  +  aso  +  2a))](aso  +  3so  +  2a)/[47V2(so  +  2)2(20so  +  3  +  0)sq], 
S  =  {3(aso  +  3so  +  2a)  (1  —  so)[(2so  +  1)(— 3sq  +  asp  —  4aso  —  12so  —  12a)0 
+18(sq  —  l)(aso  +  2sq  +  2a)]/[8-/Vs2(sp  +  2)“(20sp  +  3  +  0)2( — 1  +  4so)]}-®f, 


S2  =  Wr,  T1=- 


A1S+2B1B.+  D1T+(-2K2F1-G1K) 
Ci  > 


K  u  =  —R  —  T,  K21  =  K22  =  -K21  -S-Tu 

AT  —  Gol  ( ^21  A2  cosh(A  y)  ,  K22H 2  cosh (y,y)  .  9  q\ 

1  54 Er  y  cosh(A)  '  cosh  f. i  '  J 

/ cosh(A  y ) 


aN  sq{1—  4  so)  ex 
3 


(^r  - 1)  +  k22  - 1)  +  ^  (i/2  - 1) 

.  (-4a0-\-3so-\-16a0sQ  —  12a-\-12asQ  —  12a0so)ai  /  KnA2  cosh(Ay)  .  0  td\ 

“•  108Er  cosh(A)  h  Z 


Nso(l-\-2  so)  act 
3 


K 1 


( 


cosh(A  y ) 


coShpo  y  +R(y2~l)\ 

+  427g°a  [4a(so  -  l)(2s0  +  l)(3s0  +  1)0  +  3(-2 as0  +  3s0  +  6as02  -  4a)]y2 

+  ~a27Er  a  [^S°  +  l)(^so  —  5S0  —  4)0  —  6(s0  +  2 )\{y2  —  1)  +  2  fi2  Ks0 2  ( y 2  —  l), 

a  r  _  _  Get  (  K21  A2  cosh  A  y  ,  K22M2  cosh  fi  y  .  9  o\ 

^  5AEr  '  coshfA')  ‘  cosh  ll  '  2 


54EV 

aNso(— 1+4  so)o: 
3 


cosh(A)  1  cosh  /t, 

7^  /  cosh(Ay) 

21  - 1 


cosh(A) 

.  (8a0+12a— 3+8aso  — 16asQ  — 12so)soo:  /  K\\ 
'  1 08Er  y 


cosh  (fi) 
A2  cosh(Ay) 


)  +  *» 

2  f? 


h(A) 


2iVso(l— So)ao: 
3 


All  (3^-l)+A<!l2-l) 


+  427E°ra  [2a(l  ~  So)(2so  +  l)2)^  +  3(4aso  -  9so  -  6as02  +  2 a)y2 
+  "“27 Er  a  [—  (2so  +  l)(3«o  —  4so  —  2)0  +  3(so  +  2)](y2  —  1) 

- 2  (/xi  +  /x2)  A's02  (y2  -  1)  , 

G  =  2  ad  +  4a0Sp  +  6aso  —  3so  —  12asg  —  16a0Sg  +  lOa0so  +  6a. 
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Appendix  E.  The  coefficients  of  the  steady  solutions  for  tangential  an¬ 
choring  in  Poiseuille  flows. 


Hi  = 


6  (so +2)  (3+0(1  — so)) 

4Qtiso+377)(so+2)(3+0(l-so))+c«so(l-Ai)(-a0(l-so)(so+2)+9so(l-Az,))  ’ 

9(A_l  — 1) 


TT  _  _ 

2  ~  4(;tlso+377)(so+2)(3+e(l-so))+aso(J-Ai)(-a0(l-so)(so+2)+9so(l-Az,))  ’ 

D  =  -Ci(9iJ22E2  +  6H2F2  +  2G2HiH2), 

E  =  ~12A1(9H2E2  +  QH%F2  +  2G2HiH2)  +  C1{6H%E2  +  6H2F2  +  2G2i?i7t2) 
+12A2(9H2E1  +  6H2F1  +  2G2H±H2), 

F  =  2A1{6H2E2  +  6 HlF2  +  2G2H1H2)  -  C1H2E2  +  QH2E2 
-2  A2{QH2E1  +  6H%Fi  +  2G2H\H2), 


Rn  —  ~  §, 


T?  _  12BRn~\-E 

19  — - 


D  _  24^^11+25^12+^ 

13  — - 


12  - c - ’  ^13  - C~ 

K\\  =  (A2+2-|  [■ — 2  ^  ^ ^  —  12. -Rn  —  2 R12  +  (Rn  +  J?i2  +  +3)+], 


Kin  = 


12  (-A 2+m2) 


[4gg(A2Ai-g1.42)  _  u_Rii  _  2^  +  (jRn  +  ^  +  Rl3)X% 


ts  -  k11(b2\2+d2) 

i\oi  — - 


A2A2 


_  Ki2{B2[l2 +D2) 

-A  90  — - 


A2M2 


_  D2R12-6H2E2-6H2F2-2G2H1H2 


7?21  =  — 


Ni  = 


Ga 
54  Er 


^K2 1 


12A2 

A2  cosh(  A  y) 
cosh(A) 


_  2.B2Rl2~\-D2Rl3~\~E2H2 


,  aNso(l— 4 so)c 


7^  ^ cosh(A  y)  t  ) 
A21  ^  cosh  (A)  _  L) 


R22  =  — 

-  l2R2\y2  +  2  S22 

+  K22  ( Sfe?  - 1)  +  («“  - 1)  +  (+  - 1) 


K22/a2  cosh (y  y) 
cosh  fi 


.(—4a6-\-3so+16a6sQ  —  12a-\-12asQ  —  12a6so)<y.  (  Kn  A2  cosh(Ay)  ,  Ki2fi2  cosh (/ty) 
'  1085r  y  cosh(A)  '  cosh  (y.) 

+12Ruy2  +  2  R\ 


ATso(l+2  so)  aa 
3 


)  +A'12  (- 


cosh (y  y) 
cosh(y) 


-lj+iZn^-iJ+iZw  (y2-i) 


-P  ,/22+  [4a  (so  —  l)(2so  +  l)(3so  +  1)0  +  3( — 2aso  -P  3so  +  6aso2  —  4a)] (3j/2  —  l)2 
2asoa  #22  [(2so  + 1)  (3S2  -  5s0  -  4)0-  6(s0  +  2)]y2  (y2  -  1)  +  4M2  HxH2s<?y2  (j y 2  -  l), 


JV2  =  - 


Ga 
54  Er 


+  ypg.rf1'1'  + 12^.»2  + 2  + 


aNso(— 1+4  so)o: 
3 


All  (^?-l)+Jr”( 


cosh(y  y) 
cosh(y) 


- 1 


++i  (y4  _  i)  +  s22  (jj2  —  l) 

.  (8a0+12a— 3+8aso  — 16asg  — 12so)soQ:  /  Kn  A2  cosh(Ay)  .  Ki2/-i2  cosh  (/^y) 
'  1085r  ^  cosh(A)  '  cosh(y.) 


+  12i?ni/2  +  2  i?n 


2iVso(l— so)  act 
3 

H2soa  r 


All  (St?  -  l)  +  A12  -  l)  +  An  (+  -  1)  +  An  (+  -  1) 


1  27Er  [  2a(sQ  -  l)(2s0  +  l)(6s0  +  1)0  +  3(4as0  -  3s0  -  6os02  +  2a)](3y2  -  l)2 
^^[-(2a0  +  1)(3sq  -  4s0  -  2)0  +  3(so  +  2 )}y2{y2  -  1) 

+4  +i  +  H2)  H1s0y2  ( y 2  -  l)  , 

G  =  2  aO  +  4a0Sg  +  6aso  —  3so  —  12asg  —  16a0Sg  +  lOa0so  +  6a. 

(66) 
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Appendix  F.  The  coefficients  of  the  steady  solutions  in  normal  anchor¬ 
ing  in  Poiseuille  flows. 

tt  _  _ 6(s0+2)(9(2so  +  l)+3) _ 

1  4(/i1so+3r))(«o+2)(0(2so  +  l)+3)+a«o(l+AL)(a9(so+2)(2so  +  l)+9so(l+Ai1))> 

JLT  _ 9(1  +  Al) _ 

2  4(/iiSo+3r/)(so+2)(^(2so  +  l)+3)+aso(l+A^)(a0(so+2)(2so  +  l)+9sQ(l+A£,)) 7 

n  _  2(9HiE1+6H2F1+2G1H1H2)+(9H2E2+6HiF2+2G2H1H2) 
ni  ~  2D1+D2  7 

o  _  2(6H^E1+6H^F1+2G1H1H2)+(6H^E2+6H^F2+2G2H1H2)-12(2B1  +  B2)R1 
01  —  2Di+D2 

rji  2H2  E±-\-2H2  E2-\-2(2Bi~\-B2)Sl  jy  T  TT>  l  CT  l  'T'  \ 

Jl  = - 2D1  +  D2 - 7  Ah  =  +  *1  +  li), 


TS  _  B\A  -\-D\  is 

^21—  AiA2+Ci  n’ 


_  Difli  +  (9.H|gi+6-H^Fi+2Gig1g2) 


Ri  =  - 


(§H  %  E\+§H  %  F1+2G\H\H2)  — 12  A\R2  —  12R\B\~D\S\ 


H2  Ei+2AiS2-\-^SiBi~\-DiTi 


K22  —  ( A"21  +  i?2  +  <52  +  T2), 


_  Ga  l  K2i\2  cosh(Ay)  .  K22fi2  cosh  (jjy)  .  y2  4-  2  S  ) 

54 Er  l  cosh(A)  '  cosh '  25/  '  2  / 


cosh(A  y ) 


- 1  +  k2 


cosh (fi  y ) 


+aN,°(i,~4,j)°  h*  -  0 + **  ('sag?  -  0 

+#2  (y4  - 1)  +  s2  (y2  - 1) 

,  (— 4a0+3so  +  16a0SQ  — 12a+12asQ  — 12a0so)o:  ( Kn  A2  cosh(Ay)  ,  ir)D  rti2  ,  r>  q  A 
“•  108Er  cosh(A)  r  j 


iVso(l+2  so)  act 
3 


"tfii  (^w  ~1)+Ri  (y4  - !)  +  si  ( y 2  - 1)' 


+  g22^°a  [4a(s0  -  l)(2so  +  l)(3s0  +  1)0  +  3(— 2 as0  +  3s0  +  6as02  -  4a)](3j/2  -  l)2 
2asoa  h^sq  +  l)(3s2  -  5s0  -  4)0  -  6(s0  +  2 )}y2  (y2  - 1)  +4/z2  HiH2s02y2  ( y 2  - 1) , 

at  _  Ga  (  K2\\2  cosh  A  y  ,  K22fi2  cosh  fi  y  ,  i  o  D  ,,2  ,  o  c  A 
iV2-_54EF^ - cosh(A) -  H - ^shjl - 1H2V  +  1  *2  ) 


cosh (A  y) 
cosh(A) 


-  1  )  +  I<2 


cosh  (/j,  y) 
cosh  (/i) 


+  °A,.°14+,,")°  hi  -  l)  +  K22  -  l) 

+R2  (y4  -  l)  +  S2  ( y 2  -  l) 

,  (8a6»+12a— 3+8aso  — 16as2  — 12s0)soa  /  An  A2  cosh(Ay)  ,  10D  -.2  1  n  c  \ 

+  108Ar  ^  cosh(A)  +  i/2lT0 

+  2Nso(l-so)aa  ^  -  l)  +  R,  (/  -  l)  +  Si  ( y 2  -  l)' 

+  H|7a|°a  [— 2a(s0  -  l)(2so  +  l)(6s0  +  1)0  +  3(4as0  -  3s0  -  6as02  +  2a)](3y2  -  l)2 
~asoH2  aj_(2So  +  1)(3sq  —  4so  —  2)0  +  3(so  +  2)\y2(y~  —  1) 

+4  (m  +  /i2)  H1s0y2  (y2  -  l)  , 

G  =  2  ad  +  4a0Sg  +  6aso  —  3«o  —  12asg  —  16a0Sg  +  lOa0so  +  6a. 

(67) 
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